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ABSTRACT 

In this work we characterize the recently discovered active main belt object 
P/2012 F5 (Gibbs), which was discovered with a dust trail > 7' in length in the 
outer main belt, 7 months prior to aphelion. We use optical imaging obtained on 
UT 2012 March 27 to analyze the central condensation and the long trail. We 
find B-band and R-band apparent magnitudes of 20.96 ± 0.04 mag and 19.93 ± 
0.02 mag, respectively, which give an upper limit on the radius of the nucleus of 
2.1 km. The geometric scattering cross-section of material in the trail was ~ 4 
x 10 8 m 2 , corresponding to a mass of ~ 5 x 10 7 kg. Analysis of infrared images 
taken by the Wide-Field Infrared Survey Explorer in September 2010 reveals that 
the object was below the detection limit, suggesting that it was less active than 
it was during 2012, or possibly inactive, just 6 months after it passed through 
perihelion. We set a 1-a upper limit on its radius during this time of 2.9 km. 
P/2012 F5 (Gibbs) is dynamically stable in the outer main belt on timescales of 
~ 1 Gyr, pointing towards an asteroidal origin. We find that the morphology 
of the ejected dust is consistent with it being produced by a single event that 
occurred on UT 2011 Jul. 7 ± 20 days, possibly as the result of a collision with 
a small impactor. 

Subject headings: Comets:individual (P/2012 F5), Minor planets, asteroids: individual 
(P/2012 F5) 
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Introduction 



In recent years, a new class of objects with asteroid-like orbits but comet-like 
behavior has been identified within the main belt. These objects occupy low-eccentricity, 
low-inclination orbits with Tisserand parameters (Tj) > 3, placing them squarely in the 
asteroid regime. However, they exhibit mass loss in a manner more akin to those observed 
around comets. Though only a handful of these objects are known at this time, they cover 
a wide range of orbital space (2.29 AU < a < 3.20 AU, 0.12 < e < 0.34, 0°.2 < i < 21° .4, 
where a, e, and i represent semi-major axis, ec centricity, and inc lination, respectively), 



their sizes vary by several orders of magnitude ( IBauer et al 



20121 ). and the morphologies 



observed include spherical comae, dust tails, and/or persistent debris trails. The scientific 
community has yet to agree upon a name for such objects. For this work, we find tha t 



the terms "active asteroid" (IJewitt 



20101) , "activ e main belt object" ( 



belt comet" (IHsieh fc Jewitt 



201 2j), "activat ed asteroid" (ILicandro fc Campins 



Bauer et al. 



20121 ). and the original definition of "main 



20061 ) all apply. Potentially, depending on the driver of the 



observed activity, e ither the refined definition of main belt comet or "disrupted asteroid" 



( iHsieh et al. 



2012al ) may also apply. For simplicity, in this work we choose to utilize the 
term "active main belt object" (AMBO), which encompasses all objects with comet-like 
morphologies and main belt asteroidal orbits, regardless of the suspected driver of activity. 



Studies have shown that most o 
timescales longer than 100 Myr (e.g. 



the known AMB Os are dynamically stable on 



Hsieh et al. 



2012al lbl) . suggesting they are long-term 
residents of the main belt. Conversely, 238P and P/2008 Rl are stable for only ~ 20 
30 Mvr, resp ectivel y, and thus mavj,e interlopers from elsewhere in the Solar System 



( iHaghighipour 



2009 



Jewitt et al.ll2009h . 



One of the best characterized AMBOs, 133P/Elst-Pizzaro, has been active at multiple 
epochs, leading to the suggestion that the activity was driven by seasonal heating of an 
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(Elst et al. 


1996; 


Hsieh et al. 


2004, 


2010) 



also shown repeated activity, lending credence to the sublimation hypothesis (IHsieh et al 



20111 ). Their behavior mimics that of dynamical comets, w hich become a ctive within a 



Schorghoferl (120081 ) suggests 



few AU of the Sun as volatile deposits are heated. Work by 
that sub-surface ice deposits could survive in the main belt for billions of years at depths 
of just a few meters. Spectroscopic searches for gas emission lines have, t o date, been 



unsu c cessful but this may stem from the wea kness of the outgassing (e.g 



2009 



Licandro et al 



2011 



Hsieh et al 



Jewitt et al 



2012ah . Observations of 596 (Scheila) and P/2010 



A2 (LINEAR) did not appear to fit with the model of mass loss driven by prolonged 
sublimation. Instead, these asteroids display ed morphologies better explained by impu lsive 



mass loss events, consistent with collisions (IJewitt et al. 



2010 



Jewitt et al. 



201l|). 



Snodgrass et al. 



2010 



On UT 2012 Mar. 22.89, A. R. Gibbs reported the discovery of a new comet that 
appeared with a narrow dust trail greater than 7' in l ength. Subsequen t observations 
confirmed the discovery and the apparent morphology (IGibbs et al.ll2012l ). The object was 
given the cometary designation P/2012 F5 (Gibbs), hereafter referred to as P/2012 F5. 
On UT 2012 Jun. 21, ephemerides from the JPL HORIZONS Systeno resulted in orbital 
elements of a = 3.0038 AU, e = 0.042, and i = 9°. 739 (Tabled]), an orbital period of 
5.21 yr, and a Tisserand parameter of Tj = 3.23. The corresponding l-a uncertainties are 
6.3 x 10~ 5 AU, 8.4 x 10~ 5 , and 4°.l x 10~ 4 for a, e, and % respectively. The comet-like 
appearance but asteroid-like orbit of P/2012 F5 make it the 9 th known AMBO. 

In this work, we characterize P/2012 F5 using optical imaging data obtained at 
Palomar Observatory just 3 days after its discovery. We also use archived infrared images 
taken by the Wide-Field Infrared Survey Explorer (WISE) in September 2010 when the 



2 http:/ /ssd.jpl.nasa.gov/?horizons 
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comet was 6 months post-perihelion. Through aperture photometry, we constrain the size 
of the nucleus, and estimate the mass of dust within the trail. We use the morphology and 
dynamical modeling of the trail to characterize the duration and onset time of the activity, 
as well as the properties of the dust grains emitted. Finally, the dynamical stability of 
P/2012 F5 is investigated through the use of a symplectic integrator. Our results are 
summarized in the conclusions. 



2. Observations and Data Reduction 

In this paper, we use optical groundbased observations and infrared spacebased 
observations to characterize P/2012 F5. Table [2] provides a summary of observations. 



2.1. Large Format Camera, Palomar Observatory 



We observed P/2012 F5 with the Large Format Camera (LFC) (ISimcoe et al.l 120001 ) 
mounted on the 200" Hale telescope atop Mount Palomar on UT 2012 Mar. 27. The array 
of six 2048 x 4096 pixel CCDs provided a 24' on a side field of view, with a resolution 
of 0".35 pixel -1 when using 2x2 binning. A dithering pattern was invoked to provide 
coverage in the ~ 15" gaps between the chips. We used B and R Bessel filters with central 
wavelengths and band widths of 4400 A and 1000 A, and 6300 A and 1200 A, respectively. 
Two 180 s exposures and 5 90 s exposures were obtained with the B and R Bessel filters, 
respectively. The object was tracked at non-sidereal rates as given by the JPL Horizons 
ephemeris. 

The images were debiased and flattened using bias and twilight flatfield frames 
obtained on the night of the observations. Amplifier glow impacted all of the images in the 
southern-most chip. We investigated the extent of the contamination by comparing the 
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median of a ~ 114" x 20" box near the area of P/2012 F5's trail closest to the region of 
image most affected by the amplifier glow to the median of a background region far from 
both the amplifier and P/2012 F5. We found that the amplifier glow had no effect on the 
R-band photometry, but that the amplifier glow contributed ~ 0.03 mags to the background 
flux at the north-western end of the trail in the longer exposures obtained in B-band. 



Photometric calibration was done using Landolt standards (lLandoftlll992l ). Photometric 
uncertainties, including photon statistics, were found to be less than 0.05 mag. The debiased, 
flattened images were median-combined to produce a deeper single image in each filter. The 
stellar PSF full-width half-maximum (FWHM) was on the order of 1".7 for the individual 
images. 

The observations were obtained when P/2012 F5 was approximately 7 months 
pre-aphelion (true anomaly = 141°) at a heliocentric distance of 3.10 AU. The object 
appeared as a bright yet unresolved condensation in the south-east corner of the image with 
a long trail (Figured]). 



2.2. Wide-Field Infrared Survey Explorer 



We used archived data from the 40 cm WISE telescope to search pre-discovery images 
for signs of P/2012 F5. The data were obtained between UT 2010 Sep. 22 and 24 after 
the cryogen was depleted in the secondary tank. Consequently, only data in the 3 shortest 
wavelength bands (Wl: 3.4 yum, W2: 4.6 /im, W3: 12 jum) were available. The WISE field 
of view is 47' x 47' with a pixel scale of 2". 75 pixel -1 . Simultaneous exposures were taken 
in ea ch band every 11 s , with nominal exposu re times of 7.7 s in Wl and W2, and 1.1 s in 



W3 (IWright et al. 



2010 



Mainzer et al. 



2011ah 



Instrumental, photometric, and astrometric calibrations were performed by the "first 
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pass" scan/frame pipeline (ICutri et al.ll201l[ ). The apparent velocity of P/2012 F5 at the 
time of observation was 31" hr -1 , corresponding to a drift across the frame of less than 
0".01 during the integrations. Since this is far smaller than the pixel scale of each image, 
the effect of trailing was negligible. 

We calculated the positional uncertainty of P/2012 F5 using information retrieved from 
the JPL Horizons Ephemeris service on UT 2012 Jun. 21. The 3-cr positional uncertainty 
was 2'. 5. We thus rejected any frame in which the predicted position was less than 5' from 
the edge of the field of view. This resulted in a total of 11 useful scans of the prediscovery 
field. The individual images in each band were shifted to compensate for the motion of 
P/2012 F5 and were co-added us ing the "A WISE Astronomical Image Co-adder" (AWAIC) 



algorithm (IMasci fc Fowler 



20091 ). resulting in resampled images with pixel scales of 1".0 
pixel" 1 . The FW HM of the stacked images were 6".l, 6".4 , and 6". 5 in bands Wl, W2, and 



W3, respectively (jWright et al. 



2010 



Mainzer et al 



2011al ). The 3-cr uncertainty associated 



with the on-sky velocity of P/2012 F5 was small enough that the resulting deviation from 
its predicted position was several times smaller than the FWHM of each image. 

The observation dates of the pre-discovery fields correspond to 6 months post-perihelion 
(true anomaly ~ 38°) when P/2012 F5 was at a heliocentric distance of 2.90 AU. 



3. Results 

3.1. Nucleus Photometry at Optical Wavelengths 

Though the nucleus of P/2012 F5 is obscured by a dust coma, a bright condensation 
exists at the leading end of the trail. We use aperture photometry to constrain its brightness 
with the goal of characterizing the source of the ejecta. The aperture radius was selected 
to be ~ 2.5 times larger than the FWHM of the images but small enough to focus on 
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the material closest to the nucleus. The aperture was centered on the opto-center of the 
condensation. We note that the aperture contains non-negligible contamination from the 
coma and dust trail, and present the resulting cross-sections and radii as upper limits on 
the nucleus only. Using a circular aperture with a radius of 2".l (3260 km as projected on 
the sky at the heliocentric distance of P/2012 F5 at the time of observation), we calculate 
apparent B-band and R-band magnitudes of = 20.96 ± 0.04 mag and hir = 19.93 ± 
0.02 mag, respectively. The B-R color of the near-nucleus region is 1.03 ± 0.04 mag, which 
is consistent with solar colors. 



We correct for the observing geometry by converting the apparent magnitudes, mg, 
run, to absolute magnitudes, Hb, Hr, assuming a phase coefficient, (3 n , of 0.04 mag deg -1 , 



as is common for both active and inactive co mets (ILamv et al 



2004) as we ll as C-type 



20001 ). We calculate 



asteroids, which dominate the outer main belt (IBelskaya fc Shevchenkol 
a phase correction of 0.23 mag for the data obtained at Palomar. The corresponding 
absolute magnitudes in B-band and R-band respectively are 16.66 ± 0.04 mag and 15.63 
± 0.02 mag. The concept of absolute magnitude is not strictly valid for extended objects, 
given that an aperture of fixed angular size will include varying amounts of coma when the 



comet is observed at different heliocentric distances. 



3.2. Photometry of the Trail 

We rotate the stacked images of P/2012 F5 so that the trail is aligned horizontally 
with the image axis with the head of the trail to the left. We place a rectangular aperture 
around the trail, using the visible extent of the trail to set its boundaries. The trail is 
not symmetrical along its breadth and extends farther north-east than it does south-west. 
The photometric aperture, therefore, is set to extend 8". 75 (1.4 x 10 4 km as projected on 
the sky) in the north-east direction, and 5". 25 (8.1 x 10 3 km projected distance) to the 



-9 - 



south-west and centered on the opto-center of the trail along its length. The aperture is 
then divided into 144 (B-band) or 148 (R-band) segments along the length of the trail, 
each having dimensions of 3". 5 x 14". 0. The aperture used in B-band had a total size of 
504"x 14", while the R-band aperture had a size of 518"x 14", corresponding to physical 
distances of 7.8 x 10 5 km by 2.2 x 10 4 km, and 8.0 x 10 5 km by 2.2 x 10 4 km at the 
heliocentric distance of P/2012 F5 at the time of observation. Due to the large size of the 
trail and relatively small motion of P/2012 F5, numerous background objects contaminate 
the aperture. We digitally remove the 4 worst offenders from the stacked B-band image, 
and the 6 worst from the stacked R-band image using the IRAF task IMEDIT, which 
replaces marked pixels with values interpolated from nearby sky regions. We conservatively 
estimate that background objects remaining in the aperture contribute an uncertainty of ~ 
0.2 mag. The conversion of photometric magnitudes to physical properties, such as mass of 
the dust, are strongly dependent on the assumed properties of the dust, many of which, 
such as particle size distribution, reflectivity, and density, are uncertain by more than 20%. 
We therefore consider the intrusion of background objects an additional source of error 
that is of order or less than other sou rces. In ord e r to estimat e the sky background along 
the trail, we use a similar method to 



HsiehetaL 



(12004 



20101 ) and place rectangular sky 



background apertures of 3". 5 x 2". 8 (5.4 x 10 3 km by 4.3 x 10 3 km projected) directly 
above and below the apertures along the trail. The sky background is computed as the 
median of the pixels within these sky apertures, and is subtracted from each pixel within 
the box apertures placed along the trail. 

Figure [2] shows the normalized brightness profiles of the trail, as measured along its 
breadth and length. The cross-section through the trail's breadth is noticeably asymmetric, 
with the north-eastern edge (positive distance from the nucleus in Figure |2]) being more 
diffuse than the sharper south-western edge. We note that many cometary dust tails and 
trails show asymmetries, perhaps most noticeably the trail of P/2010 A2 (LINEAR), which 
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retained an unusual and highl y asymmetric morp 



an ou tburst in the main belt (IJewitt et al. 



2010 



lology for 9 month s after undergo-in 



Moreno et al 



2010 



g 



Snodgrass et al 



20101 ) . Thus, this asymmetry may be a result of the intrinsic structure present in the 
ejected material shortly after the outburst. The surface brightness decreases along the 
length of the trail with increasing distance from the nucleus, d, approximately as d~ 0A . 
This is a somewhat shallower decrease t han the relat i onship of d~ os that was found for 



133P/Elst-Pizarro as observed in 2002 by 



Hsieh et al. 



(120041 ) . The material within the trail, 



as measured using the rectangular apertures previously described, has total magnitudes 
of 17.41 ± 0.20 mag and 16.25 ± 0.20 mag in B- and R-band, respectively. We use 3". 5 
x 14". segments to measure the B - R color along the trail between the nucleo-centric 
distances of 14" and 504". We find the median color of the segments to be 1.12 ± 0.28 
mag, consistent with both the color of the region at the head of the trail, and solar colors. 
We note that the large uncertainty al so renders it consistent with the color of dust around 



active comets (B-R~1.0-1.4 mag; 



Jewitt fc Meech 



1986 



Kolokolova et al 



200J) 



3.3. Constraints From WISE Infrared Photometry 

Using the stacked WISE images, an area of 600" x 600" was searched for any PSF-like 
signal that exceeded a l-a limit above the sky background using a 9".0 radius aperture (1.8 
x 10 4 km projected distance). P/2012 F5 was not detected in any wavelength band by the 
WISE mission. By characterizing the sky background in W3 (the most sensitive band to 
thermal dust emission) we set a 3-cr detection limit on any potential objects within the field 
of 5 mJy, which corresponds to a lower limit on the absolute magnitude of material within 
the aperture of 15.2 mag. 
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4. Analysis 



4.1. Dust Morphology 



The standard method for unders tanding the morpho 



trails is the Finson-Probstein model ( IFinson fc Probstein 



ogy_ 



of comet ary dust tails and 



1968). This model assumes that 



once cometary dust particles leave the surface, their motion is governed by two forces: solar 
gravity and solar radiation pressure. The particle motion can then be parameterized using 
the ratio, (3 of these two forces: 



In physical units, this gives the ratio: 



rad 



grav 



1.19 x l(r 3 Q 



pr 



Pdd 



(2) 



where Q pr is the scattering efficiency for radiation pressure, pd is the bulk density of 



the particle, d is the particle diamet er, and the factor o 



multiplyi ng all the constan t values (IFinson fc Probstein 



1.19 x 10 3 [kg m 2 ] comes from 



Q pr ~ 1 (IBurns et al. 



1968). For grain sizes of d > A, 



19791 ): given the central wavelengths of the filters used, this condition 



holds for d > 0.5 pm. 



j3 is incorporated into the equation of motion in the following way: 



, „, GMp, _ 
x +(1-/3)— = 
\xr 



(3) 



where G is the universal gravitational constant, M & is the mass of the Sun, x is 
the vector position of the object. This is a simple equation of motion that can then be 
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integrated for different values of f3 to track the motion of particles with a particular /3 value. 



Th e computations w ere carried out by creating a numerical integrator (based on the 



work of 



LisseetaL 



19981 ) in the language Python which took in a set of /3 values (0.0001 
< (3 < 3.0000), and integrated the motion of the dust particles over the designated time 
interval. This generated a set of points which can be shown as curves of particles with 
constant j3 released at a range of times (syndynes) or curves of constant release date with a 
range of particle sizes (synchrones). 

Figure [3] shows plots of the syndyne models plotted on top of the R-band data. At first 
look, the (3 = 0.03 (shown in magenta) syndyne seems to model the trail (highlighted with 
a white dashed line) well. However, upon closer inspection it is clear that this syndyne 
initially is north-east of the trail, then curves south-west and crosses the trail approximately 
200" (projected distance: 3.1 x 10 5 km) from the nucleus. 

When the model is instead shown as synchrones (curves of constant particle release 
date, Figure Hj), the fit is greatly improved. The particle emission date was constrained by 
modeling the width of the trail with a Gaussian, which yielded a half-width half-maximum 
of about 1".05. We then calculated where the /3 = 0.03 syndyne came within this distance 
from, and then crossed, the center of the trail. This revealed a best-fit synchrone of 264 ± 
20 days, as shown as a green line in Figure HI This corresponds to a particle ejection date 
of UT 2011 Jul. 7 ± 20 days when P/2012 F5 was at a heliocentric distance of 3.01 AU and 
a true anomaly of ~ 94°. 

We use the observed length of the trail to constrain the size of the particles, given that 
the distances travelled by dust grains since the ejection event are size- dependent. We find 
that particles with radii of ~ 20 /im could have traveled ~ 8 x 10 5 km to the end of the 
observed trail in 264 days, setting a lower limit on the size of particles present. Smaller 
particles could have traveled further but are not observed in our data, possibly because 
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they were either not originally released or had a scattering cross-section that fell below our 
detection limit. 

4.2. Physical Characterization of P/2012 F5 and its Activity 

The results from the aperture photometry can be used to constrain the physical 
properties of the nucleus and adjacent coma. We use Equation H] to estimate the geometric 
cross-section of material, a [m 2 ] within the aperture. 



$„a = 2.24 x 10 22 vr r^A 2 10 a4(m ^- m ^ 



PX^a 



(4) 



where p\ is the geometric albedo, m Q \ is the apparent magnitude of the Sun, and rri\ 
is the apparent magnitude of the material within the aperture, each given for a specific 



broadband filter ( iJewitt 



19911 ) . We assume the albedo of material is 0.04 in both B- 



and R-bands, which is broa dly consistent with observations of AMBOs, outer main be 



asteroids, and comets (e.g. 



Masiero et al. 



2011 



Jewitt & Meech 



Bauer et al. 



1986 



Lamy et al 



2004 : 



Hsieh et al. 



2009a 



t 



20121 ). The apparent magnitudes of the Sun in B-band and 
R-band are -26.10 and -27.12, respectively Q. The parameter $ a is a function to correct for 
phase-angle-dependent variations in brightness and is calculated using Equation [5j 



- 2.5log 10 <& a = a/3 



(5) 



where a [deg] is the phase angle and (5 a [mag deg l ] is the phase coefficient, already 
defined as having a value of 0.04 mag deg -1 . 



1 http://mips. as. arizona. edu/^cnaw/ sun. html 
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Considering the apparent magnitudes of material within a 2".l aperture, we estimate 
the geometric cross-section of the material to be ~ 1.4 x 10 7 m 2 in B-band and ~ 1.5 x 10 7 
m 2 in R-band. If the light was being reflected by an inactive, spherical body, the effective 
radius of the object would be ~ 2.1 km. This sets an upper limit on the size of P/2012 F5, 
though, given the extensive nature of the coma, the radius of the asteroid is probably on 
the order or less than 1 km. 

We also use the non-detection of P/2012 F5 in the WISE data to set an upper 
limit on the size of a bare nucleus that could have been present during September 2010. 



The upper limit on the flux from the W3 sta cked image and the NEATM model ([Harris 



998 



Delbo et al. 



( jLebofsky et al. 



2003 



Mainzer et al. 



2011bl ) with a fixed beaming (rj) parameter of 0.8 



19861 ) and an assumed optical albedo of 0.04 ± 0.02 yielded a radius of < 



2.9 km, for a 1-cr confidence interval. 

The upper limits on the size of P/2012 F5 are consistent with reported sizes of other 
AMBOs, which range in radius from 120 m to 113 km, with a median radius of 2.2 km. It 
should be noted that all but one AMBOs have reported radii under 5 km. 

To determine whether P/2012 F5 would have been observed by WISE if it displayed 
a similar brightness and morphology to its appearance in the Palomar images, we compare 
the lower limit on absolute magnitude, H R WIS e ^ 15.2, to that of P/2012 F5 during 
the Palomar observations. We use a 11". 9 radius aperture (projected distance 1.8 x 10 4 
km) to match the physical distance subtended by the 9".0 radius aperture used with the 
WISE data. Centering this aperture on the condensation observed in March 2012, we 
find an apparent R-band magnitude of 18.67 ± 0.02 mag, corresponding to an absolute 
magnitude of 14.37 ± 0.02 mag. This is ~ 0.8 mag brighter than the limit set by WISE, 
suggesting that P/2012 F5 would have been detectable by WISE in September 2010, had it 
displayed 2012-levels of activity. We can thus say that the object must have had a different 
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morphology when observed in 2010, and was possibly a bare asteroid. 



We use the qua ntity Afp as a measure o 



A'Hearn et al. 



(119841) and further discussed in 



formulation given in 



A'Hearn et al. 



dust product i on as initially defined in 



A'Hearn et al. 



(119951 ) . We adapt the initial 



(119841 ) to give: 



Afp = p (^J^j wO-H^-m) (6) 

where A is the geometric albedo, / is a filling factor of the grains in the field of view, 
p is the projected size of the aperture [cm], A is the geocentric distance [cm], r H is the 
heliocentric distance [AU]. m is the phase-corrected magnitude of the material, while m Q x 
retains its previous definition. The model assumes the coma follows simple radial outflow 
model, which is insufficient to explain the observed morphology of P/2012 F5. We thus 
restrict ourselves to using photometric results from small apertures (< 3") when calculating 
Afp. Since the coma does not follow the assumed radial model where N(p) oc p, where N(p) 
is the number of dust grains within an aperture of size p, we find that the calculated values 
of Afp are not independent of aperture size. We find that the value of Afp peaks at p ~ 
the FWHM of the image, and subsequently decreases by ~ 1 cm for every 1" increase in 
aperture radius. For an aperture of radius 2".l (a projected distance of 3260 km), Afp = 
11.0 cm and 11.1 cm in B-band and R-band, respectively. The uncertainties are dominated 
by the morphological divergence of the data from the model of simple, symmetric outflow 
and are significant enough that the results from the different filters are consistent with each 
other. This suggests that the grains are gray in color. 

Equation H] can also be used to estimate the geometric cross-section of material within 
the trail if the apparent magnitudes of the trail are substituted for m\. We find a = 3.6 x 
10 8 m 2 and 4.1 x 10 8 m 2 for the B- and R-band magnitudes, respectively. As discussed in 
Section 14.11 the minimum particle size within the trail is ~ 20 pm, while the largest dust 
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grains close to the nucleus may be > cm-sized. By making the simplifying assumption that 
the bulk of the mass of the trail resides in spherical grains with radii of 100 /im and bulk 
densities of 1000 kg m -3 , we find that the mass of material within our apertures is ~ 5 x 
10 7 kg. Taking the radius of P/2012 F5 to be ~ 1 km with a bulk density of 1000 kg m -3 , 
the material in the trail amounts to ~ 0.01% of the asteroid by mass. 



4.3. Dynamical Stability 



We next seek to constrain the origin of P/2012 F5 to test whether it it likely to have 
originated in the main belt, or to have been recently inserted from elsewhere in the Solar 
System. Using numerical simulations to model test particle evolution, we have investigated 
the stability of 100 clones of P/2012 F5 on the Gyr-timescale. Clone starting positions were 
randomly assigned according to a Gaussian distribution with a FWHM equal to the 1-cr 
uncertainties in the osculating elements of P/2012 F5 given by the JPL HORIZONS System. 
We integra ted the orbits of each clone for 1 G yr using the SWIFT_RMSVY symplectic 



integrator (ILevison fc Duncan 



1994 



Broz 



20061 ) with a 25 day timestep, and including the 



Venus, Earth, Mars, Jupiter, Saturn, Uranus, and Neptune as massive bodies in addition 
to the Sun. We ignore non-gravitational effects for the purpose of this simulation as they 
are under-constrained. No clones were lost from the system during the simulation and 
the orbital evolution of the clones is negligible over significant periods of time (Figure E]), 
indicating that P/2012 F5 is dynamically stable and is not likely to have been recently 
implanted in the Main Belt. 
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5. Discussion 



Dynamical modeling of the ejecta of P/2012 F5 suggests that the activity was driven by 
a one-time event, approximately 260 days before our observations. This rules out continuous 
sublimation over a period of months, though a short, intense burst of sublimating material 
over a period of days could recreate the morphology observed. Considering the possibility 
that the activity is driven by the sublimation of sub-surface volatiles, we examine the 
orbital location of P/2012 F5 during the WISE and Palomar observations. P/2012 F5 
was at a heliocentric distance of 2.90 AU and 6-months post-perihelion during the WISE 
observations, which revealed that the object was not displaying the activity observed 
at its discovery in early 2012. Strangely, P/2012 F5 was discovered only months before 
reaching aphelion, at a time when its surface temperature is decreasing as its heliocentric 
distance increases. However, the eccentricity of P/2012 F5 is only 0.04, corresponding to 
only a ~ 6 K difference in surface temperature between perihelion and aphelion (assuming 
it absorbs and emits as an isothermal blackbody with an albedo of 0.04). We therefore 
do not expect P/2012 F5 to be active only near perihelion, if the activity is driven by a 
temperature-dependent process, such as sublimation. It is difficult to reconcile the short 
period of activity with sublimation when dust tails, trails and/or comae observed around 
other AMBOs, such as 133P, 238P, and P/ 2010 R2, are cons i stent w i th prolonged generati on 



over a period of several months or more (IHsieh et al. 



2004 



2009b 



2012b 



Moreno et al. 



2011al ). We note that if the estimated ejected mass lost through sublimation over a period 
of two weeks, the mass loss rate would be ~ 25 kg s _1 , orders of magnitude higher than 
that observed for other AMBOs whose activity is suspe cted to be driven by sublimation 



(IHsieh et al. 



2004; 



Jewitt et al. 



2009 



Hsieh et al. 



2009bJ). Observations of repeated activity 



on P/2012 F5 would strengthen the argument for sublimation-driven activity, especially if 
subsequent mass loss occurred at a similar orbital position to the 2011 outburst studied 
here. 
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Smaller astero ids are subject 



rotational-spin up ( IRubincam 



;o non-gravitational forces that can cause rapid 
2000). If the rotational velocity of m aterial on the surface 



exceeds the escape velocity, mater ial may be ejected 



not well constrained, but work by 



Walsh et al. 



((Jewitt 



2012l ) . The fate of ejecta is 



(120081 ) suggests that the material would be 



lifted into a low orbit where it may eventually accrete into a satellite. This is supported 
by observations of asteroid 1999 KW4, which is rotating with a period close to the limit, 



(Ostro et al. 


2006 


Scheeres et al. 


2006) 



20061 ). Due to a lack of well-studied examples, it is unclear 
whether rotational breakup could recreate the morphology observed around P/2012 F5, 
or whether it would eject sufficient mass on a short timescale of just a few weeks. Future 
measurements of the rotation period and size of P/2012 F5 could assist with gauging the 
likelihood of this scenario by calculating whether the nucleus is rotating rapidly enough to 
eject material. 

The results of the dust modeling are also consistent with an impact that occurred ~ 
260 days prior to our observations. If the activity observed around P/2012 F5 is due to a 
collision with a smaller asteroid, the mass of observed ejecta can be used to set an upper 
limit on the size of a possible impactor. A mass of 5 x 10 7 kg, as calculated in Section I4.2[ 
would correspond t o a sphere of radius 20 m, assuming a bulk density of 1000 kg m~ 3 . As 



Jewitt et al. 



(120101 ) note, the majority of the mass lost from a small body during an impact 
is from the primary body, rather than the impactor. We thus conclude that the radius of 
an impactor would have been on the order of meters in size. The ej ecta mass is comparab le 



to estimates of the ejected mas s around P/201 A2 (3.7 x 



5 x 10 7 kg, 



loLkg, 



Moreno et al. 



Hsieheta 



201 ll : 6 x 10 8 kg, 



2012a 



20 id 6-60 x 10 7 kg, 



4 x 10 7 kg, 



Jewitt et al. 



Jewitt et al 



2011 



Bodewits et al. 



2011 



2 x 10 iU kg 



7 kg , ISnodgrass et al. 



2010 



2010) and 596 



Scheila) (3 x 



1.5-4.9 x 10 s kg 



Moreno et al 



Ishiguro et al 



201 lbl ). both of which 



are suspected to have undergone collisions in the main belt. With so few examples, it is 
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difficult to know if the similarity of the ejecta masses is coincidence or intrinsic to the 
nature of main belt asteroidal collisions. 



6. Summary 

We have used optical and infrared data to characterize P/2012 F5, which appears to 
be a main belt asteroid undergoing mass loss as a result of an impact with a smaller body. 

1. A 1-a upper limit on the radius of the nucleus is set at 2.9 km using WISE observations 
from 2010. A more sensitive upper limit of 2.1 km is set using R-band images from 
Palomar observatory. The cross-section of material observed near the nucleus in 2012 
is above the detection limit set by the WISE data, suggesting that the asteroid was 
either inactive or weakly active during September 2010. 

2. Activity ejected ~ 3 x 10 8 m 2 of material into the trail, spanning > 7' (6.5 x 10 5 km 
projected distance) in space. The minimum size of particles in the trail was ~ 20 ^m. 

3. We find that syndynes do not fit the observed morphology, arguing against a period 
of continuous mass loss from the nucleus. Instead, a single event on UT 2011 Jul. 7 ± 
20 days can explain the observed shape of the trail. 

4. Results from a symplectic integrator show the object to be dynamically stable and 
resident in the outer main belt on long timescales, suggesting P/2012 F5 is native to 
the main belt. 

5. The observed behavior can be explained by a collision with a meter-sized impactor. 
Further observations may reveal morphological evolution in the trail that may shed 
more light on the cause of the activity. 
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Table I. Orbital Elements of P/2012 F5 



Epoch 


a [AU] 1 


q [AU] 2 e 3 


i [deg] 4 


fi [deg] 5 


w [deg] 6 


M [deg] 7 


2012 Mar. 28.0 


3.00 


2.88 0.04 


9.74 


216.9 


177.4 


138.4 



1 Semi- major Axis 

2 Perihelion 

3 Eccentricity 

4 Inclination 

5 ongitude of the ascending node 

6 Argument of perihelion 

7 Mean anomaly 



Table 2. Journal of Observations 



Date (UT) Observatory Filter N ! Exposure time [s] r H [AU] 2 A [AU] 3 a [deg] 4 Image scale [km pixel x ] 



2010 Sep. 22-24 


WISE 


W3 


11 


1.1 


2.90 


2.7 


20.2 


5400 


2012 Mar. 27.3 


Palomar 


R 


5 


90 


3.10 


2.1 


5.8 


535 


2012 Mar. 27.3 


Palomar 


B 


3 


180 


3.10 


2.1 


5.8 


535 



1 Number of exposures 

2 Heliocentric distance 

3 Geocentric distance 
4 Phase angle 
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Fig. 1. — P/2012 F5 as observed on UT 2012 Mar. 27 using the Large Format Camera 
mounted on the Hale 200" telescope atop Mount Palomar. The image was created from a 
median stack of 5 90 s R-band exposures taken while tracking the motion of P/2012 F5. The 
velocity vector, v, and the solar direction, 0, are marked, as are the directions North and 
East. 
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Fig. 2. — R-band surface brightness profiles of P/2012 F5 along the breadth (left) and 
length (right) of the trail, calculated using a ~ 500" x 14" (8 x 10 5 km x 2 x 10 4 km 
projected) aperture. Uncertainties due to photon statistics are plotted, though the scatter 
is increased through contributions from background objects. Especially bright background 
stars are labeled with 'S'. 
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Fig. 3. — Syndynes tracing particles with (3 = 0.0003 - 3.0 (bottom to top, cyan through 
yellow) are plotted over a stacked R-band image of P/2012 F5 from UT 2012 Mar. 27. The 
syndynes do not recreate the observed morphology and diverge from the trail (highlighted by 
the white dashed line) at large nucleo-centric distances. The best-fitting syndyne (/3 = 0.03, 
magenta) initially curves north of the trail then crosses to the south-west approximately 200" 
(3.1 x 10 5 km projected) from the nucleus. We find that synchrones provide a much better 
fit, suggesting the dust was released during a single event, rather than over a prolonged 
period. 



Fig. 4. — The best-fit synchrone (green) is plotted on the stacked R-band image with addi- 
tional synchrones representing the uncertainties (red). The position and approximate length 
of the trail is marked by the dashed white line. Particles represented by the best-fit syn- 
chrone must have been released 264 ± 20 days prior to the observation, which corresponds 
to an ejection date near UT 2011 Jul. 7. 



T = Gyr 



T = 1 Gyr 




i>+ + + 
+ t + + + 




3.0028 3.0033 3.0038 3.0043 3.0048 3.0028 3.0033 3.0038 3.0043 3.0048 2.9988 3.0013 3.0038 3.0063 3.0088 2.9988 3.0013 3.0038 3.0063 3.0088 
a[AU] a [All] a [All] a[AU] 



Fig. 5. — Dynamical evolution of 100 clones (black crosses) that started with orbits within 
the quoted 1 a uncertainties of the orbit of P/2012 F5 (filled red circle). After precessing 
the orbit for 1 Gyr, the particles have dispersed in orbital space but not significantly. The 
object P/2012 F5 can thus be considered a dynamically stable main belt asteroid, and not 
an interloper from the outer Solar System. 



